ON EIGENVALUES OF LAME OPERATOR 



KOUICHI TAKEMURA 

Abstract. We introduce two integral representations of monodromy on Lame 
equation. By applying them, we obtain results on hyperelliptic-to-elliptic reduc- 
tion integral formulae, finite-gap potential and eigenvalues of Lame operator. 



1. Introduction 
The Lame operator is a differential operator defined by 

where is the Weierstrass p-function and n is a constant. In this manuscript we 
assume that basic periods of the function p{x) are (1, r) and is a positive integer. 

We are going to consider eigenvalues of the Lame operator H with boundary con- 
dition. In other words, we specify the vector space where the operator H acts and 
consider the spectral problem on the space. 

Let f{x) be an eigenfunction of the operator H with eigenvalue i.e. 



H = -—+n{n + l)p{x). (LI) 



{H - E)f{x) =(^-^+n{n + l)p{x) - f{x) = 0. (L2) 

Then Eq. ()1.2|) is called the Lame equation (or the Lame's differential equation). 
Although the differential equation p.2|) is periodic in x with period 1, the solution may 
not be periodic. We also note that the eigenfunction may not be square-integrable 
on the interval (0,1) because the potential n(n -|- l)p{x) has poles on Z. In this 
manuscript we consider the condition for the eigenvalue E such that Eq. ()1.2|) has a 
singly-periodic eigenfunction (i.e. /(x-|-l) = ±/(x)), a doubly-periodic eigenfunction 
(i.e. f{x+l) = ±f{x) and f{x+r) = ±/(x)) or a square-integrable eigenfunction (i.e. 
lo \fi^)\'^d,^ < +oo). These conditions correspond to the boundary condition and the 
boundary condition is closely related with the monodromy. For the case n G Z, the 
monodromy of solutions to Eq. p.2|) is calculated and it has two expressions. One 
is expressed by hyperelliptic integral and the other is based on Hermite-Krichever 
Ansatz. Note that hyperelliptic-to-elliptic reduction integral formulae are obtained 
by comparing two expressions. 

A crucial fact on the Lame equation is that the Lame operator is an example of 
the finite-gap potential which was established by Ince 0. By applying the formula of 
the monodromy in terms of hyperelliptic integral, we can show results on finite-gap 
potential. Moreover relationship between the boundary condition and the finite-gap 
potential is clarified. 
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The eigenvalue E of the Lame operator with boundary condition depends on the 
period r. In particular we can numerically investigate branching of the eigenvalues 
as a complex-analytic function in r for the case n = 1, and it is compatible with the 
convergence radius of the eigenvalue expanded as a power series in p = exp{n\/—lT). 

Here we comment on relationship with quantum integrable system. The elliptic 
Calogero-Moser-Sutherland model (or elliptic Olshanetsky-Perelomov model ^T]) of 
type An is a quantum many body system whose Hamiltonian is given as follows, 

1 92 

^^iv = "2 5Z + ^^'^ + ^) J2 Pi^i-^j)- (1-3) 

i=l * l<i<j<N+l 

This model is known to be integrable, i.e. there exists iV+l-algebraically independent 
commuting operators which commute with the Hamiltonian H^^. For the case = 1, 
we recover the Lame operator by setting x = Xi — X2 and y = Xi + X2- Therefore the 
elliptic Calogero-Moser-Sutherland model is a generalization of the Lame operator. In 
contract to the case of the trigonometric Calogero-Moser-Sutherland model, spectral 
problem of the elliptic Calogero-Moser-Sutherland model is not studied very much. 
We hope that the spectral problem of the elliptic Calogero-Moser-Sutherland model 
is clarified by developing knowledge on the Lame equation in future. The paper by 
Chalykh, Etingof and Oblomkov |2j might suggest us an approach to this problem. 

This manuscript is organized as follows. In section El we describe relationship 
among Lame equation, Heun's differential equation ^2] and hypergeometric polyno- 
mial. In section 121 we introduce a doubly-periodic function that satisfies a differential 
equation of order three, and obtain an integral representation of solutions to the 
Lame equation and a formula of global monodromy in terms of hyperelliptic integral. 
In section HI we explain results on Bethe Ansatz and Hermite-Krichever Ansatz. As 
an application of Hermite-Krichever Ansatz, we get a formula of global monodromy 
in terms of elliptic integral. In section we obtain hyperelliptic-to-elliptic reduc- 
tion integral formulae by comparing two expressions of monodromy. In section (HI 
we describe results on the finite-gap potential. In section [71 we consider analytic 
continuation of eigenvalues in the variable r for the case n = 1. 



2. Relationship with other equations 

It is known that Lame equation is a special case of Heun equation. The Heun 
equation (or the Heun's differential equation) is a differential equation given by 

+ ^ + "fi^-^ ,, )/W-o (2.1) 

aw J \w w — 1 w — t J aw w[w — l)[w — t) I 
with the condition 

j + S + e = a + (3 + l. (2.2) 

Note that Eq. ()2.1|) has four singularities {0,l,t, 00} on the Riemann sphere and 
the Heun equation is the standard canonical form of a Fuchsian equation with four 
singularities. It is well known that the Fuchsian equation with three singularities is 
the hypergeometric differential equation. 
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The Heun equation admits an expression by elhptic function. Let Hjj be the 
Hamiltonian of BCi Inozemtsev model which is given as 

(f ^ 

i^iT = -^ + 5^/i(/i + l)p(a; + cu,), (2.3) 

i=0 

where uq = 0, ui = 1/2, U2 = — (r + l)/2 and = r/2 are half-periods, and k 
{i = 0,1,2,3) are coupling constants. This model is a one-particle version of the 
BCjsf Inozemtsev system [Hj, which is known to be the universal quantum integrable 
system with B]\r symmetry 01101 • Let /(x) be an eigenfuction of the operator 
with eigenvalue E, i.e. 

Hnfix) = Ef{x). (2.4) 

Under the transformation 

w = — , w—iw-l) — {w-t) — f{w) = f{x), (2.5) 

p(x) - 63 

Eq. ()2.4p is transformed to the Heun equation ()2.H) . In this sense, Eq. ()2.4|) is an elliptic 
representation of the Heun equation. The coupling constants Iq, li, I2, I3 correspond 
to the exponents a, . . . , e, the elliptic modurus r corresponds to the singular point t 
and the eigenvalue E corresponds to the accessory parameter q. For details see [T3] . 
If /o 7^ and h = h = h = 0, then Eq. ()2.4p represents the Lame equation. Thus the 
Lame equation is a special case of the Heun equation. 

We observe a relationship with hypergeometric polynomial. More precisely, hy- 
pergeomertic differential equation is obtained from the Lame equation (or the Heun 
equation) by trigonometric limit. Set 



HT = -^ + n{n + l)-^. (2.6) 

ax^ sm TTX 



Then H — > — '^n{n + 1) by the trigonometric limit J9 = exp(7rv— Ir) 



3 

The equation {Ht — E)f{x) = is transformed to the hypergeometric differential 
equation. Set 

<|)(2;) = (sin7rx)"+\ = CmC^i^^ {cos ttx)^{x), (m G Z>o), (2.7) 
n^,F,(-m,m + 2.;. + i;i^ 



where the function C^(z) = ^^^^2^) 2Fi{—m,m + 2iy; + |; ^-g^) is the Gegenbauer 
polynomial of degree m and = ^J '^ '''r(^+2n+2) ""^^ • Then we have 

HTVm = 7r^{m + n + iyVm, (2.8) 

and the set {vm}m&i^o forms a complete ortho normal system on a certain Hilbert 
space. Hence the spectral problem for the operator Ht is clarified substantially. 

To investigate the spectral problem for the Lame operator H, we can apply a 
method of perturbation from the operator Ht and we have an algorithm for obtaining 
eigenvalues and eigenfunctions of the operator H as formal power series in p. It is 
shown in ^Hl HI that, if \p\ is sufficiently small, then the formal power series 
of eigenvalues and eigenfunctions converge. Relationship between the convergence 
radius and the monodromy will be mentioned later. 
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3. MONODROMY AND HYPERELLIPTIC INTEGRAL 

In this section we review results on an integral representation of solutions and a 
monodromy formula which is expressed by hyperelliptic integral. For this purpose we 
introduce a doubly-periodic function which plays important roles. Let h{x) be the 
product of any pair of the solutions to Eq. (jl.2j) . Then the function h{x) satisfies the 
following third-order differential equation: 

(P d \ 

— - 4 (n(n + l)p{x) -E) — - 2n{n + l)p'(x) h{x) = 0. (3.1) 
dx-^ dx J 

It is known that Eq. p.ljl has a nonzero doubly-periodic solution for all E. 

Proposition 3.1. [j T^ Proposition 3.5] Ij n G Z>i, then Eq. \3. i)) has a nonzero even 
doubly-periodic solution S(x, -E), which has the expansion 

n-l 

E{x,E) = co{E) +Y,bj{E)p{xr-\ (3.2) 

j=0 

where the coefficients Cq{E) and bj{E) are polynomials in E, they do not have common 
divisors and the polynomial Cq{E) is monic. 

For the case of Lame equation, we have deg^ cq{E) = n and deg^; bj{E) = j. 
We can derive an integral formula for the solution to Eq. ()1.2|) by using the doubly- 
periodic function E{x,E). Set 

Q{E) = Eix, Ef {E - n{n + l)p{x)) + ls(x, E)^^^^ - \ (^^^) ' ■ 

(3.3) 

Then the r.h.s. is independent of x, and Q{E) is a monic polynomial in E of de- 
gree 2n + 1. The following proposition on the integral representation of solutions is 
obtained in P^: 

Proposition 3.2. ^14, Proposition 3.7] Let E{x,E) be the doubly-periodic function 
defined in Proposition 1,9.11 and Q{E) be the monic polynomial defined in Eq. ( I,?.,*^) . 
Then the function 



A(x,E) = v/HO^exp/ v^^^' 
is a solution to the differential equation M.^) . 



(3.4) 



Note that Eq. ()3.4|) is written in the book of Whittaker and Watson [20j for the 
case of the Lame equation. 

Example 1. (i) The case n = 1. 

S(x, E) = p{x) + E, Q{E) = {E + ei){E + e2){E + eg). (3.5) 
(ii) The case n = 2. 

E{x, E) = 9p(a;)2 + ?,Ep{x) + E^ - \g^, Q{E) = {E^ - Sg^) UliiE - 3e,), (3.6) 
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where g2 = -4(eie2 + 6263 + 6361). 
(in) The case n = 3. 

E{x, E) = 22hp{xf + AhEp{xf + 6(^2 - '-^g^Mx) + - Ihg^E - ^g^, (3.7) 
Q{E) = E nli(^' + Qe,E + 45e2 - 15<72), (3.8) 
where g^ = 4616263. 

Now we consider the monodromy. It follows from Eq. ()3.4|) that A{x + 1, E) (resp. 
A{x + T,E)) is expressed as miA{x,E) (resp. mT-A{x, E)) for some constant mi 
(resp. rrir). We determine the constants mi and rrir. Rewrite the function S(a;,i?) 
as follows: 

H(x,E) = c(E) + ^a,(E) - pix). (3.9) 
i=o ^ 

Proposition 3.3. |Tni Assume n G Z>i. Le^ £" 5e the eigenvalue satisfying 
Q{E') = 0. Then there exists qi, qs G {0,1} such that A{x+2uJk, E') = {-iy^=A{x,E'), 
and we have 

A(x + 1,E) = {-irAix, E) exp [ T "^^^ ~ rfE | , (3.10) 

V \J-QiE) I 

A{x + T,E) = (-l)-A(x, E) exp [ T | , (3.11) 

where rji = ({1/2), rj^ = C{t/2) and C,{x) is the Weierstrass zeta function. 

Example 2. (i) The case n = 1. Set E' = —61. Then we have gi = 0, gs = 1. The 

polynomials c{E) and aQ{E) are calculated as 

c{E) = E, ao{E) = 1. (3.12) 
(a) The case n = 2. Set E' = \fig2- Then qi = qz = 0, 

c{E) = E^- 1^72, ao(^) = 3E. (3.13) 
(Hi) The case n = 3. Set E' = 0. Then qi = qs = 0, 

c{E)=E^-fg,E-'-fgs, ao{E) = Q{E^ - g,) . (3.14) 

4. Bethe Ansatz and Hermite-Krichever Ansatz 

We review validity of Bethe Ansatz and Hermite-Krichever Ansatz for the Lame 
equation, and provide examples of the Hermite-Krichever Ansatz. 

The method "Bethe Ansatz" appears frequently in physics. In this manuscript, we 
use "Bethe Ansatz" in a somewhat restricted sense. We assume that the eigenfucntion 
of the Lame operator has a special form as Eq. ()4.1|) . Then the Bethe Ansatz method 
replaces the problem of finding eigenstates and eigenvalues of the Hamiltonian with 
a problem of solving transcendental equations for a finite number of variables which 
are called the Bethe Ansatz equations (see Eq. ()4.3|) ). We review a proposition that 
almost all eigenfunction is expressed as the form of Bethe Ansatz. For the case of 
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Lame equation it was discussed in [201, extended to the case of the Heun equation 
in |Sl E] • Note that it is apphed to show vahdity of the Hermite-Krichever Ansatz. 

Proposition 4.1. (i) For the case Q{E) ^ 0, there exists ti, . . .in and A such that 
tj 7^ tj' (j 7^ j'), tj ^ + |Z and A(x, E) is expressed as 

A(x, E) = A \-^^n exp -x V ({tA , (4.1) 

where o"(x) is the Weierstrass sigma function, 
(a) The function 

n? 1 <^ix + tj) 

m = -^r^ ^exp(cx), (4.2) 

(^i^r Uj=i (^{tj) 

with the condition tj ^ tj/ (j ^ j') and tj ^ |Z + |Z is an eigenfunction of the Lame 
operator (see Eq. M.l]) ). if and only iftj (j = 1, . . . ,n) and c satisfy the relations, 

- + -tk)=c (j = 1, . . . ,n). (4.3) 

The eigenvalue E is given by 

n 

E = -c'-n J2iP(tj) - C(tjf) + J2^Pitj - tk) - - tkf). (4.4) 

i=l j<k 

Now we deal with the Hermite-Krichever Ansatz. In our situation, the Hermite- 
Krichever Ansatz asserts that the differential equation has solutions that are expressed 
as a finite series in the derivatives of an elliptic Baker- Akhiezer function, multiplied 
by an exponential function. More precisely, we are going to find solutions to Eq. ()1.2|) 
of the form 

'n— 1 




fix) = exp (kx) \yj:i{^] Hx, a) ] (4.5) 

where $(x, a) = exp{({a)x)a{a — a;)/(cr(x)cr(a)). From Eq. ()4.5p we have 

f{x + l,E) = exp(-2r/ia + C(a) + k), (4.6) 
fix + T,E)= exp(-2r/3a + r(C(a) + n)). (4.7) 

If the value p{a) is calculated, then the values a and C{a) are expressed as elliptic 
integrals. Hence, if a solution /(x) is expressed as Eq. (|4.5p . and p{a) and k are 
calculated explicitly, then we obtain a monodromy formula by elliptic integrals. 

Hermite and Halphen investigated the Lame equation for the cases n = 1, 2, 3, 4, 5 
on this Ansatz in the 19th century. Belokolos, Eilbeck, Enolskii, Kostov and Smirnov 
studied for the case Iq = 2,li = 1, Z2 = 0, Z3 = and the case Iq = 2,li = 1, 12 = 
l,/3 = on the Heun equation (see and the reference therein). On the other 
hand, Treibich and Verdier jTHj constructed a theory of elliptic soliton following 
Krichever's idea and studied "tangential covering" which would be closely related 
with the Hermite-Krichever Ansatz. 

Now we remind validity of the Hermite-Krichever Ansatz. 
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Theorem 4.2. ^H] Assume n G Z>i. There exists polynomials Pi{E) , . . . , P^^E) 
such that, if P2{E) ^ 0, then a solution to Eg. il.^) is written as Eg. Ii4.5\ ) for some 
values a, k, and bj (j = 0, . . . , n — 1). The values a and k are expressed as 

Note that a in Eq. ()4.5|) and tj in Eq. ()4.H) satisfy the relation a = — Yl^=i ^ji which 
is crucial to show validity of the Hermite-Krichever Ansatz, and the function A(x, E) 
in Eq. ()3.4|) is also expressed as the r.h.s. of Eq. ()4.5|) . 

To calculate the polynomials Pi{E) , . . . , P^i^E) , it is effective to apply notions 
"twisted Lame polynomial" and "theta-twisted Lame polynomial" which were in- 
vented by Maier They are generalized to the case of the Heun equation in [TH] . 

Example 3. (i) The case n = 1. 

p{a) = -E, K = 0. (4.9) 

(a) The case n = 2. 

{E -3ei){E + 6eif 2 



^(")=^^ 9iE^-3g2) - "= 3(E^-3,.) ^-^("^^- ^''''^ 



Note that the term {E + 6ei) corresponds to the twisted Lame polynomial 
(Hi) The case n = 3. 

(^2 + 6eiE + Abej - 15^2)(^^ + ISciE - 2256^ + f c/2) 
SQE{E^ - f (?2 



= ei o^^,^, 7^^^ , (4.11) 



K= ^^^^Q{E). (4.12) 

6E(E2-f(?2)^ ^ ^ ^ ^ 

Note that the terms {E'^ + 15eiE — 22'be\ + ^fi's) (ind {E"^ — ^(72) correspond to the 
twisted Lame polynomials. 

5. HYPERELLIPTIC-TO-ELLIPTIC REDUCTION INTEGRAL FORMULA 

We obtained two expressions of the monodromy in Eqs. ()3.1()| IH.llj) and Eqs. ()4.(i^ 
14. 7|) . By comparing these two expressions we obtain hyperelliptic-to-elliptic reduction 
integral formulae. For the proof, see 



Theorem 5.1. [IS., §3] Let Pi{E), . . . .P^i^E) he the polynomials defined in Theorem 
\4-S\ and let ao{E),c{E) be the ones in section{^ Set 

(i) We have a formula 



2 



dE= / ^ =, (5.2) 



which reduces a hyperelliptic integral of the first kind to an elliptic integral of the first 
kind. 
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(ii) Let E' be the value that satisfies Q{E') = and P2{E') ^ 0, and let ^' be the 
value ^ (see Eg. 115.1]) ) evaluated at E = E' . Then we have a formula 



1 r c{E) 



-.dE = -K + 



-Q{E) 



^d^ 



(5.3) 



4^ - 92^ - 93 



which reduces a hyperelliptic integral of the second kind to an elliptic integral of the 
second kind. 



Example 4. (i) The case n = 1. The values ^ and n are given by ^ 

Hence the formulae 45.4 trivial. 

(ii) The case n = 2. The values ^ and n are given by 



-E and n = 0. 



(E-3ei)(E + 6ei 



9{E^-3g2) ' 3{E^-3g2)\ 



iE^-3g2)lliE~3ei) 



i=l 



(5.4) 



Set E' = 3ei. Then we have ^' = ci. Eqs. \5.^A are written as 



3E 



-{E^-^g^)Y[^^^^{E-?,e, 



-.dE 



^2 3 _ 
E - 2^2 



3ei 



-.dE = -K + 



^ \/-(^'-3^2)n-=i(^-3e.) 



4^' - 92^ - 93 
' \/ 4:^-92^ -93 



(5.5) 
(5.6) 



These formulae reduce hyperelliptic integrals of genus two to elliptic integrals. By 
setting ^ = —y/Q, E = z, g2 = a/3, g^ = 6/54, we recover the formula 



zdz 



dy 



^/{z'^-a){8z^-6az-b) 2^3 J ^y^-3ay + b' \ 3{z^ - a 



2z^ - b 



(5.7) 



from Eg. / I5.,5|) . This formula was found by Hermite in the 19th century. From Eg. h5. 6) 
we have 



{2z'^ — a)dz 



1 Sz'^ — Qaz — b 1 



ydy 



y/{z^ -a){8z^ -6az-b) SV - a 2^3 J y/y^-3ay + b 

(Hi) The case n = 3. The values ^ and k are given by 



. (5.8) 



(^2 + QeiE + 45e? - 15^2) (^^ + 15ei^ - 225e? + f 



3QE{E^ - f (72)2 



K 



QE{E^ - f ^2) \ 



-E Y[{E^ + GciE + 45ef - 15^2). 



(5.9) 
(5.10) 



i=l 
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Let E' be a root of an equation + GciE + 45e^ — 15(72 = 0. Then we have ^' = Ci. 
Eqs. \5.'A \5.c^) are written as 

-\ /\ . ^.^^^-^.^-^ t , (5.11) 



-^nti(^' + 6e.i? + 45e2- 15(72) sj - 92^ - 9, 



^' V-^nti(^' + 6e.i? + 45e2 - 15^2) ^^^^ ^J - 92^ - ^73 

(5.12) 

These formulae reduce hyperelliptic integrals of genus three to elliptic integrals. 

6. Boundary value problems and Finite-gap potential 

We consider boundary value problems of the Lame operator H. Let aint {H) be the 
set of eigenvalues of H whose eigenvector is square-integrable on the interval (0, 1), 
i.e. 

a,nt{H) = {E I 3f{x) ^ s.t. Hf{x) = Ef{x), [ \f{x)\dx < +00}. (6.1) 



Let <7ci{H) be the set of eigenvalues of H whose eigenvector is doubly-periodic, i.e. 

a,{H) = {E I 3/(x) ^ s.t. Hf{x) = Ef{x), f{x + 1) = ±/(x), f{x + r) = ±/(x)}, 

(6.2) 

Note that the doubly-periodic eigenvector is simply the Lame polynomial up to gauge 
transformation and changing variable. It is known that ^ad{H) = 2n -|- 1 and 

a^iH) = {E I Q{E) = 0}. (6.3) 

Let (Ts{H) be the set of eigenvalues of H whose eigenvector is singly-periodic, i.e. 

a,{H) = {E I 3/(x) ^ s.t. Hf{x) = Ef{x), f{x + 1) = ±f{x)}, (6.4) 
On the sets aint{H), (yd{H) and crs{H), we have 
Proposition 6.1. [17j Assume that r G V~lI^>o- Then 

aint{H)U^diH)=as{H). (6.5) 
Next, we briefly explain a relationship with finite-gap potential. Let 

Hs = -^ + n{n + l)p{x + T/2) (6.6) 

be the operator which is obtained by the shift x —>■ x + t/2. Then the potential of 
the operator Hg does not have poles on R and it is real-analytic, if r G V~lI^>o- Let 
ab{Hs) be the set such that 

E G (TbiHs) -v^ All solutions to {Hg — E)f{x) = are bounded on x G M. 

It is known in the theory of Hill's equation jH| that, if the potential of the op- 
erator —d^/dx^ + q{x) is real-analytic and singly-periodic, then the open set M \ 
<yb{—d'^/dx'^ + q{x)) is expressed as (—00, Eq) U Ui(zi{E2i-i, E2i) where / = {1, . . . , n} 
{n G Z>o) or / = Z. If the number of the unstable bands (gaps) is finite (i.e. 
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I = {1, . . . ,n} for some (n G Z>o)), then the potential q{x) is called finite-gap po- 
tential or finite-band potential. 

It is known that the potentials of the Lame operator (see Eq. fl6.6|) ) for n G Z>i is 
typical examples of the finite-gap potential, which was established by Ince |5j. 

Theorem 6.2. Assume that r G V— lM>o and n G Z>i. Then 



M \ ab{Hs) = (-00, Eq) U (El, E2) U ■ ■ ■ U (E2„_i, ^sn) (6.7) 

where Eq < Ei < ■ ■ ■ < E2n and Ei G <yd{H) (i = 0, . . . , 2n). 

Now we present a brief sketch of a new proof of the theorem which is based on 
the monodromy formule in terms of hyperelliptic integral. From Eq. ()3.10|) we have 
A(x + 1, E) = ± exp(— mi/2)A(a;, E), where 

-Q{E) 



E' 



and E' satisfies Q{E') = 0. Note that, if r G v^K>o, then ao{E),c{E),Q{E) G M 
for G M. If Q{E) > 0, then the integrand is pure imaginary by choosing E' 
appropriately, and we have an equality | exp(— mi/2)| = 1. Thus the eigenfunction is 
bounded. If Q{E) < 0, then the integrand is real by choosing E' appropriately, and 
unboundedness of the eigenfunction is shown. 

In the theory of stationary soliton equation (see ^ and the reference therein), alter- 
native definition of the finite-gap potential is known. If there exists an odd-order dif- 
ferential operator A = [d/dxf'^^^ +X]^lo^^i(^) {d/ dx)^^~^~'^ such that [A, —d'^/dx^ + 
q{x)] = 0, then q{x) is called the algebro-geometric finite-gap potential. Note that the 
equation [A, —d'^/dx'^ + q{x)] = is equivalent to the function q{x) being a solution to 
some stationary higher-order KdV equation. It is known that, if q{x) is real-analytic 
on M and q{x + 1) = q{x), then q{x) is a finite-gap potential if and only if q{x) is an 
algebro-geometric finite-gap potential. For the case of the Lame equation (and the 
Heun equation), the operator A is calculated by using the function E{x,E). Write 

n 

H(x,E) = ^a„_,(x)E\ (6.9) 
From Proposition 13. II we have do{x) = 1. 

Proposition 6.3. UHI Define the {2n + l)st-order differential operator A by 

A = ± [a,(. + r/2) A _ ^ (^a,(. + r/2)) } ff-. (6.10) 

Then the operator A commutes with the operator Hg. Moreover we have 

A^ + Q{H,) = Q. (6.11) 
Another formula for the operator A of determinant type is obtained in ^21 §3.2]. 
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7. The case n = 1 and analytic continuation of eigenvalues 
In this section we consider the case n = 1. The Hamiltonian is given as 

^ = -^ + 2K^)- (7-1) 

Eigenvalues which admits doubly-periodic eigenfunction are written as 

ad{H) = {-ei,-e2,-es}. (7.2) 

The functions Ci, 62 and 63 are analytic on the upper half plane M + a/— l]R>o a^^d 
they are related with automorphic functions. 

Next we consider eigenvalues which admits square-integrable eigenf unctions. Set 
p = exp(7r-\/— Ir). Then r V— loo corresponds to p 0, and r G M + ■\/~lI^>o 
corresponds to \p\ < 1. For the case p = 0, a complete orthonormal basis for square- 
integrable eigenfunctions are written as ijn G Z>o) in Eq . (12.71) with the eigenvalue 
Em = vr^(— 2/3 + (m + 2)^). Based on the eigenvalues Em and the eigenfunctions 
Vm for the case p = 0, we determine eigenvalues Em{p) = Em + YlT=i ■^"^''^ P^'^ 
normalized eigenfunctions Vm{p) = f m + Zlfcli Sm'ez>o ^m%'"rn'P'^^ for the operator 
H as formal power series in p by applying a method of perturbation. For details see 
[TH| ITT] . Convergence of the formal power series of eigenvalues in the variable p was 
shown in [TKj . 

Proposition 7.1. [15, Corollary 3.7] Let Emip) {fn G Z>o) (resp. Vmip)) be the 
formal eigenvalue (resp. eigenfunction) of the Hamiltonian H. If \p\ is sufficiently 
small then the power series Em{p) converges. 

Note that it was also shown in [15] that the eigenvalues Em{p) are real-analytic in 
p2G(-l,l). 

The inferred convergence radius and expansions of the first few terms for the eigen- 
value Emip) are calculated as follows (see also [TT): 



Eo{p) 




' li) 1 80 2 1 1360 4 1 20800 6 , 195920 8 , 3174880 10 , \ 
^o~ro/^ ~r 27 P 243 ^ 2187 ^ 19683 ^ i ■ ■ ■ J 


.749 


E2{P) 




Ub 1 272^2 1 198928^4 | 55403584^6 | 4307155408^^,8 | \ 


.749 


E,{p) 


vr^ 


^ll)6 1 %2 2 1 i?i]248 4 , 1?'?1^33184 6 ■7^M'^62851088 8 , \ 
^ 3 3.^ 4287.5 P ' 5252187.5 P ' 128678.59375 P "T • • • j 


.875 


Eiip) 


TT^ 


+ 20f + 65/ + / + -'^I'p' + %^'p'' + -'f,tp'' + ...) 


.838 


Esip) 




'fi 1 52 2 1 1493 4 1 35671 6 , 4492153 8 , 55853449 10 , \ 
^ ' "i P ' 27 P 486 P 34992 P 629856 P ' • • • ) 


.838 


E,{p) 




'•m 1 82 2 1 50339 4 1 13640101 6 1 3872868499 's \ 
k 3 ^ f,P ^ \mnP 1 2onono P ' .32ononon P ^ ■■■) 


.906 



Table 1. Expansion and inferred convergence radius. 



If the function Eq{p) is analytic on the upper half plane M + V— 1M>05 then the 
radius of convergence in p should be 1, thought it is numerically inferred to be .749. 
Hence it seems that the function Eo{p) has singularity. 

Now we consider analytic continuation of the eigenvalues Em{p) [m G Z>o) in the 
variable p. From Eq. ()3.10j) or Eqs. ()4.6U4.9|) . the eigenvalue E is continued analytically 
in p by keeping the conditions 



E = -pito), 

3m G Z, — C(^o) + 2r7ito = mn\^^ 



(7.3) 
(7.4) 
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Let 3fJa (resp. Qa) be the real part (resp. the imaginary part) of the number a, and 
let Ca be the cycle starting from 3fJa, approaching the point a parallel to the imaginary 
axis, turning anti-clockwise around a and returning to 3?a as shown in Figure 2. 



Im 


r 

Ca 


\ 




Ra 



Re 



Figure 2. Cycle Ca- 

We continue the eigenvalue E analytically along the cycle Ca on where a is a possible 
branching point for \a\ < .90, 3fJa > and Q'a > 0. The branching points for \a\ < .90, 
5Ra > and S^a > are calculated in jT2j, and they are hsted in Table 3. Branching 
along the cycle Ca is determined as shown in Table 3 (see also [T7j). 



a = .258666 + .697448/ 
a = .224582 + .842777/ 
a = .552288 + .677536/ 
a = .314813 + .821858/ 
a = .686317 + .559106/ 


Eq{p) E2{p), E2{p) Eo{p), Ei{p) Ei{p), Ee{p) E^{p) 
Eo{p) Ei{p), E2{p) E2{p), E4p) Eo{p), E^ip) Ee{p) 
Eo{p) => E^ip), E2{p) E2{p), E^ip) => Eo{p), Ee{p) => Ee{p) 
Eo{p) => E^ip), E2{p) E2{p), E^ip) => Eo{p), Ee{p) => Ee{p) 
Eo{p) => Eo{p), E2{p) Ei{p), E^ip) => E2{p), Ee{p) => Ee{p) 


a = .535905 + .640487/ 


Ei{p) Esip), E3{p) Ei{p), E^ip) E^ip), Er{p) E^ip) 



Table 3. Branching along the cycle Ca 

The closest branching point from the origin is p = .258666 + .697448/ {\p\ = 
.743869) and the eigenvalues Eq{p) and /?2(p) are connected by continuing analytically 
along the cycle Cp {p = .258666 + .697448/) (see Table 3). On the other hand, by 
the method of perturbation the convergence radii of the expansions of the eigenvalues 
Eo{p) and E2{p) are both inferred to be around .749 from Table 1. Thus, convergence 
radii calculated by different methods are very close and compatibility between the 
method of perturbation and the method of monodromy is confirmed. Moreover, we 
obtain a reason why the convergence radii of the eigenvalues Eo{p) and E2{p) are 
very close by considering the branching point. For details see |17j . 

It would be able to consider analytic continuation of eigenvalues for the case n > 2 
and the case of the Heun equation by applying results on the Hermite-Krichever 
Ansatz. 
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